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INTRODUCTION 
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In  a  recent  paper  [1]  it  was  shown  that  there  are  smooth,  nonlinear,  three- 
dimensional  controllers,  not  incorporating  probing  signals,  which  are  capable  of 
adaptively  stabilizing  any  single-input,  single-output,  minimum  phase,  relative 
degree  two  or  less  linear  system  of  any  dimension.  Controllers  of  this  type  are 
based  on  minimal  dynamic  compensator  synthesis  [2],  While  such  controllers  are 
simple  in  structure  they  do  not  have  a  model-following  capability. 

r~-f 

In  this  paper  we  develop  a  new  algorithm  based  on  observer  theory  1-2],  which  can 
adaptively  stabilize  and  achieve  model-following  as  well.  The  controller,  which 
is  a  smooth  nonlinear  dynamical  system  of  dimension  4(n+l),  can  adaptively  sta¬ 
bilize  any  physical  process  with  scalar  input  u  and  scalar  output  y  ,  provided 
the  process  can  be  modelled  by  a  strictly-proper,  minimum  phase,  linear  system  of 
dimension  not  exceeding  n  and  relative  degree  not  exceeding  two.  The  controller 
is  based  on  concepts  developed  previously  in  [1]  and  [3],  /j 

1.  CONTROL  EQUATIONS 

The  controller  to  be  examined  consists  of  a  two-dimensional  reference  system 
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where  A^  and  A  are  positive  constants  and  r ( * )  is  a  bounded,  differentiable 
reference  input,  a  tracking  error 


e  *  y  “  yr»  U) 

sensitivity  function  n-vectors  6^  and  6^  generated  by  the  equations 
8  »  A8  +  bu 

u  u  (3) 

8  -  A8  +  by 

y  y 

where  n  is  a  prespecified  positive  integer  and  (A,b)  is  any  n-dimensional  control¬ 
lable  pair  with  A  stable,  and  a  control  law 


u  “  N(||kJ|)(k^8u  +  k’8  +  kpp  +  kfr) 


where  k^.k^.k  and  kf  are  control  parameters, 

k  »  [k* ,k',k  ,k  ] ' , 
u  y  p  r 


JJkJ|*  (k'k)  ,  and  N(*)  is  a  Nussbaum  Gain,  i.e.  any  integrable  function  satisfying 
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{e.g.,  N(p)  -  p  cos  (m  )  •  Control  parameters  k  ,k  ,k  and  k  are  adjusted 
according  to  the  rules  u  y  p  r 
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The  controller  defined  by  ( 1 ) — ( 7 )  may  be  viewed  as  a  smooth  4 (n+1 )-dimensional 

dynamical  system  with  inputs  r,  r  and  y,  state  {y  ,p,e  ,6  ,z  ,z  ,z  ,z  }  and  output 
u.  ruyuypr 

Remarks; 

1.  By  using  minimal  dimensional  observer  theory,  it  is  possible  to  reduce  the 
dimensions  of  8^  and  6^  to  (n-1)  and  to  eliminate  one  control  parameter  thereby 

obtaining  a  (4n+l)-dimensional  algorithm  with  the  same  capabilities  as  the  one 
described  here  [2],  The  stability  analysis  of  the  lower-dimensional  algorithm 
is  essentially  the  same  as  the  analysis  which  follows. 

2.  There  are  two  different  ways  to  avoid  generating  the  derivative  of  the  refer¬ 
ence  signal  r  required  by  the  above  algorithm.  The  first  is  simply  to  intro¬ 
duce  a  new  reference  signal  r  and  then  make  r  a  state  of  a  three-dimension¬ 
al  reference  system  defined  by  (1)  and  r  +  X^r  -  r  where  X^  >  0.  The  second 

is  to  make  use  of  an  idea  due  to  Monopoli  [4,5]  which  amounts  to  replacing  (4) 
and  (7)  by  ' 

u  ■  N(||k||)k’0  +  +  Nq|kj|H'*)e 


k.  —  q»e 

respectively  where 


V  3  ' 


6  -  [e^.p.r]'  (8) 

and  p  +  «  0.  This  alternative,  however,  requires  N ( • )  to  be  dif  ferentiable.}r 

2.  MAIN  RESULT  41 


The  process  to  which  the  preceding  algorithm  is  applicable  must  admit  a  minimum 
phase,  linear  model  of  dimension  n  £  n  and  relative  degree  £  2.  It  is  known  [6] 
that  these  assumptions  imply  that  such  a  process  can  also  be  modelled  by  a  n- 
dlmensional,  stabilizable,  minimum  phase,  relative  degree  two  or  one  linear  system 
of  the  form 


(A  +  h  c)x  +  b  u 
P  P  P 
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where  c  is  any  row  vector  chosen  so  that  (c,A)  is  observable,  and  k  and  b  are 

P  P 

unknown  constant,  parameter  vectors.  We  take  (9)  as  the  model  of  the  process  to 
be  controlled.  Our  main  result  is  as  follows. 

Theorem  lj  For  each  initial  state  and  each  differentiable  input  r,  bounded  on 
[Ot'Ot  the  state  response  X  m  (x^,8^,6^ ,p,y^ ,z^,z^ ,z^ ,z  )  of  the  closed-loop  system 

defined  by  (l)-(9)  exists  and  Is  bounded  on  [0,°>)  and  the  tracking  error  e  -*■  0  as 
t  -*  «.  ~  -  ~ 

The  remainder  of  this  paper  is  devoted  to  a  proof  of  this  theorem. 

3.  STABILITY  ANALYSIS 

To  prove  Theorem  1,  it  proves  useful  to  work  with  a  certain  system  of  equations 
which  we  now  derive.  For  this  let  n*  denote  the  relative  degree  of  (9)  and 
define 


6  -  y  +  Xxy 

It  is  known  [6]  that  for  n*  -  1, 


(10) 


6  -  gu  -  d  '6  -  d'0  +  e  (11) 

uuyy 

where  g  is  a  nonzero  constant  -  the  "high-frequency  gain"  -  d^  and  d^  are  unknown 

parameter  vectors,  and  c  is  an  unknown  linear  combination  of  decaying  exponentials. 
Similarly,  for  n*  ■  2  it  is  known  that 


6  +  X,6  -  gu  -  d’8  -d’0  +  e  (12) 

2  °  u  u  y  y 

where  g,  d^.d^  and  e  have  the  same  interpretations  as  for  n*  *  I.  It  is  also 
known  [6]  in  either  case  that  if  x  *  [8^.8,,x,],,  then  x  satisfies 

x  -  Ax  +  b(gu  -  d'0  -  d'0  +  e)  (13) 

uuyy 

where  A  is  stable.  * 


To  proceed,  define  % 

o  -  e  +  Xxe  (14) 

and  note  front  (1),  (2)  and  (10)  that  o  “  6  -  p.  It  follows  from  (1),  (11)  and  (12) 
that 

o"gu-d'e  -d'8  +e-p 
°  uuyy 

if  n*  •  1  and 

o  +  X0o  *  gu  -  d'0  -  d'0  +  e  -  r 

2  u  u  y  y 

if  n*  »  2.  Thus  if  we  define  d  -  [d'.d'.l.O]’  for  n*  ■  1  or  d  -  [d',d’,0,l]’  for 
n*  -  2,  then  u  y  u  y 


u  ■  gu  -  d’0  +  e  if  n*  ■  1  (15) 

and 

o  +  X2O  ■  gu  -  d'8  +  e  if  n*  ■  2  (16) 

where  0  was  defined  previously  in  (8).  If  in  addition  we  define  q  «  [1,0]'  for 
n*  *  1  or  q  *  [0,1] '  for  n*  ■  2,  then  for  n*  *  1  or  2  (13)  becomes 

x  ■  Ax  +  b(gu  -  d*8  +  c  +  (p,r]q) 


(17) 


Thus  with  the  notation 


u»gu-d'0+e, 

(15)  through  (17)  can  be  rewritten  as 


if  n*  -  1 


a  +  X^o  *■  u  if  n*  ■  2 


x  -  Ax  +  b(u  +  [p,r]'q)  (17)' 

By  differentiating  the  expressions  for  k^.k^.k^  and  kf  in  (7)  and  substituting  in 

the  expressions  for  z^,z^,z^,z^  also  in  (7),  it  is  straightforward  to  verify  that 

k  -  6  a,  k  ■  0  a,  k  *  po  and  k  -  ra,  a  being  given  by  (14).  Using  (5)  and 
u  u  y  y  p  r 

(8)  we  can  thus  write 

k  ■  0o  (18) 

For  ease  of  reference  we  now  summarize  in  one  place,  the  system  of  equations  to  be 
analyzed. 


e  +  X^e  ■  o 


a  +  X^a  -  u 


if  n*  ■  1 


if  n*  -  2 


u  ■  gN ( J |  kj|)k'e  -  d'9  +  e  (19d) 

k  -  6o  (19e) 

• 

x  -  Ax  +  b(u  +  [p,r]'q)  (19f) 

P  +  X2p  -  r  *  (19g) 

0  -  [O^.Oj.P.r]'  (19h) 

x  -  [0',0',xT  (19i) 

u  y  p 

The  preceeding  defines  a  dynamical  system  of  the  form  Z  “  F(Z,r,e)  where 
Z  -  [e,o,k,x,p]#  (with  o  deleted  from  Z  if  n*  -  1).  Observe  that  boundedness  of  Z 
implies  boundedness  of  e,o,k,0u,0^fXp  and  p.  Boundedness  of  p  together  with  (1) 

implies  boundedness  of  y  .  In  addition,  boundedness  of  e,k,0  ,0  ,p  and  r  together 

r  y  y 

with  (5)  and  (7)  imply  boundedness  of  z^jz^.z^  and  z^.  Thus  to  prove  Theorem  1 
it  is  enough  to  show  that  Z(t)  exists  and  is  bounded  on  [0,»)  and  that  e(t)  +  0  as 


Differentiability  of  F  implies  that  for  any  initial  state  Z(0),  there  must  exist 
an  Interval  I  m  (0,t^)  on  which  a  solution  Z(t)  exists.  For  any  function  5(t) 


defined  on  I,  we  write  S  c  (“(l)  if  (,  ?tis  bounded  on  l^by  a 

depending  on  t^;  we  also  write  {  c  1  (I)  if  j  ||((r)||  dr  e  1  (I). 


a  constant  not 


Vie  can  now  state 

Proposition  1:  There  exist  constants  C,  and  C„,  not  depending  on  t..  such  that  for 


J  0(T)3(T)dT  *  g*(  ||k(t)  ||)  +  J  a2(T)dr  +  k(t)jj  +  C2  (20) 


*(C)  -  pN(p)dp  (21) 

2  0 

Proof:  Set  n  “Ilk  II  ;  hence  from  (19e),  n  -  2k'8o.  It  follows  from  this,  (19d) 
and  (19e)  that  *' 

ou  «  |  N ( n*5) n  -  d'k  +  oe  (22) 

We  can  now  develop  bounds  for  the  Integrals  of  the  terms  on  the  right  side  of  (22). 
First  note  that  .  . 

ft  fn(t) 

g/2  j  N(nJi(T))n(T)dr  -  g/2  j  N(nSdn 
0  n(0) 

Ikon 

-  g  uN(ui)du 

||mo)|  (2) 

Next  observe  that  “  gtr ( ||  k( C )  J  J  )  -  gs(  ||k(0)  ||  ) 

-[  d'k(x)dx  *  d*k(0)  -  d'k(t) 


ft 

-  d'k 

Jn 


||d  I  k(t)||+  |d'k(0)  | 


2  2 

For  any  positive  constant  C,  o  e  jCo  +  e  /4C.  Thus 
[  o(x)e(x)dx  i  C  o2(x)dx  +  HC  f  e2(t)dT 

Jo  Jo  Jo 

By  setting  C  ■  A2/(l+A2),  Cj  >||d||and 

C2  -  I«'(||k(0)||)|  +  (d’k  (0)  |  +  >tc  j  e2(x)dx, 
(23)-(25)  can  be  combined  to  yield  (20).  V 
Proof  of  Theorem  1:  From  (19b)  and  (19c)  It  follows  that 


fc  2  ft2 

a  (x)dx  If  n*  ■  1  or  A  o  (x)dx 

■*0  1  Jo 

Proposition  1  thus  imply  that 

("ufe r,  f° 
«||Mt,||)  2  J  1  ,0 

/  f‘ 
L  <i+v  Jo" 


t 

o(x)u 
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(x)dx  equals 


or  A  o2(x)dx  +  *i(o2(t)  -  o2(0))  it  n*  “  2.  This  and 

“  r\ 


if  n*  ■  1 


2 

2 ( x )dx  +  y-(t)  if  n*  -  2 


0(0  -  g*(C)  +  c^.  +  C2  +  0  (0)/2 

In  view  of  (6)  and  the  definition  of  n(C)  in  (21)  it  is  easy  to  see  that  there  must 
exist  a  closed-bounded  interval  (a,bj  containing  ||k(0)[|  for  which  both  f:(a)  and 


fl(b)  are  negative.  Since  (26)  implies  that  for  any  t^e  I,  n(  j|k(t)|J)  5  0,  ||k(t)IJ 
cannot  pass  through  either  a  or  b.  Therefore JJ kjj e,L“( I) .  In  addition,  since  1 
£1(0  is  continuous,  it  follows  from  (26)  that  c  e  l  (I)  for  n*  -  1  and 

o  e  L2(I)ftL*(I)  for  n*  -  2.  Thus  a  e  12(I)  for  n*  -  1,2  so  by  (19a) 

e  e  12(I)Ql“(I)  for  n*  -  1,2. 

For  n*  ■  1,  u  «  o  e  i2(I);  since  [p,r]'  e  L  (I)  it  follows  from  (19f)  that 
x  e  L°°(I). 

For  n*  -  2,  u  »  a  +  X20  and  0  e  ^  (I),  again  it  follows  from  (19f)  that 
x  e  1*(I). 

At  this  point  we  have  shown  that  Z  -  [e,o,k,x,p]'  e  L  (I)  for  n*  «  2,  that 
'  •  80  2 

I  Z  •  [e.k.x^p]1  e  i  (I)  for  n*  *  1  and  that  e  e  L  (I]  f^r  n*  -  1,2,  Therefore  we 

can  take  t^  “  “.  Thus  Z  is  bounded  on  [0,“)  and  eel  [0,“).  In  addition,  since 

(19)  implies  that  e  e  L  (0,<“),  it  follows  that  e  -*■  0  as  t  ■*  “.  7 

CONCLUDING  REMARKS 

>  The  algorithm  presented  here  and  its  subsequent  stability  analysis  rely  for  the 

*  most  part  on  ideas  developed  previously  in  [1],  [3]  and  [6].  In  fact  the  stability 
analysis  given  is  almost  identical  to  that  used  in  [1].  The  essential  new  idea  in 
this  paper  is  to  use  a  control  law  (4)  incorporating  both  reference  model  state  p 
and  reference  input  r.  It  is  this  departure  from  more  traditional  adaptive  control 
laws,  (e.g.  [3])  which  makes  model  fallowing  possible  with  one  algorithm  for 
processes  of  both  relative  degree  1  and  2. 
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